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Article  history:  Ocean  thermal  structure  is  critical  for  the  intensity  change  of  tropical  cyclones  (TCs).  It  has  been  opera- 

Available  online  26  October  2013  tionally  derived  from  satellite  altimetry  for  TC  forecasting  and  research.  The  existing  derivation  is,  how¬ 

ever,  based  on  a  simple  two-layer  method;  as  a  result,  only  two  isotherms  can  be  obtained  to  coarsely 
characterize  the  subsurface  ocean  thermal  structure.  Improvement  on  the  vertical  resolution  to  enhance 
ocean  characterization  is  desirable  for  more  accurately  estimating  ocean’s  energy  supply  for  TC  intensity 
change. 

In  this  study,  we  present  a  new  generation  of  derivation  to  improve  ocean’s  subsurface  characterization 
for  the  Western  North  Pacific  Ocean  (WNPO)  because  this  region  has  the  highest  TC  occurrence  on  the 
Earth.  In  addition  to  the  presently  used  two  isotherms  for  the  depths  of  20  °C  and  26  °C  isotherms 
(D20  and  D26),  we  derive  continuous  isotherms  from  D4  up  to  D29  (maximum  26  subsurface  layers) 
to  characterize  the  subsurface  ocean  thermal  structure  in  detail.  This  is  achieved  through  applying  a  large 
set  (>38,000)  of  in  situ  Argo  thermal  profiles  to  regression  development.  A  smaller  set  of  in  situ  Argo  pro¬ 
files  (>7000),  independent  of  those  used  for  regression,  is  utilized  for  validation,  to  assess  the  accuracy  of 
the  new  derivation.  The  root-mean-square  differences  (RMSDs)  between  the  derived  and  the  in  situ  iso¬ 
therms  are  found  to  be  within  ~10-20  m  for  the  upper  isotherms  (D20  to  D29)  and  within  ~40-60  m  for 
the  lower  isotherms  (D4  to  D19).  No  significant  biases  of  derived  isotherms  are  found. 

In  addition  to  increasing  the  vertical  resolution  from  two  layers  to  26  layers,  higher  accuracy  for  the 
two  isotherms  of  D20  and  D26  is  also  achieved.  For  example,  in  the  existing  two-layer  method,  D20  in 
the  northern  WNPO  is  grossly  overestimated  with  a  high  RMSD  of  ~80-100  m;  the  new  method  signif¬ 
icantly  reduces  the  RMSD  to  ~20  m.  Better  subsurface  characterization  leads  to  improvement  in  impor¬ 
tant  TC-intensity  related  parameters,  like  during-cyclone  sea  surface  temperature  (SST)  cooling  and  air- 
sea  enthalpy  flux  supply.  Based  on  a  series  of  ocean  mixed  layer  numerical  experiments  using  16  ran¬ 
domly-selected  profiles,  we  find  that  the  error  in  SST  cooling  (air-sea  flux  supply)  can  be  reduced  from 
36%  (120%)  to  13%  (20%). 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Ocean  thermal  structure  (OTS)  plays  a  critical  role  in  the  inten¬ 
sity  change  of  tropical  cyclones  (TCs),  due  to  the  energy  supply 
from  the  ocean  to  the  TC  (Leipper  and  Volgenau,  1972;  Gray, 
1979;  Holliday  and  Thompson,  1979;  Emanuel,  1999;  Schade  and 
Emanuel,  1999;  Shay  et  al.,  2000;  Goni  and  Trinanes,  2003;  Lin 
et  al.,  2005,  2008,  2013;  Wu  et  al.,  2007;  Goni  et  al.,  2009;  Shay 
and  Brewster,  2010;  Balaguru  et  al.,  2012;  Moon  and  Kwon, 
2012;  Clone  et  al.,  2013).  From  the  perspective  of  TC,  OTS  typically 


*  Corresponding  author  at:  Department  of  Atmospheric  Sciences,  National 
Taiwan  University,  No.  1,  Sec.  4,  Roosevelt  Rd.,  Taipei  106,  Taiwan.  Tel.:  +886  2 
33663917. 

E-mail  address:  iilin@as.ntu.edu.tw  (l.-l.  Lin). 

0079-6611/$  -  see  front  matter  ©  2013  Elsevier  Ltd.  All  rights  reserved. 
http://dx.d0i.0rg/l  0.1 01 6/j.pocean.201 3.1 0.004 


concerns  temperatures  in  the  top  200-300  m,  because  the  interac¬ 
tion  between  TC  and  the  ocean  takes  place  within  this  layer.  When 
a  TC  passes  over  the  ocean,  the  underlying  OTS  is  disturbed  and 
mixed  by  the  intense  wind  associated  with  the  TC.  As  a  result, 
the  relatively  cold  water  from  deeper  ocean  is  entrained  and/or  up- 
welled  into  the  surface  layer,  producing  the  renowned  TC-induced 
sea  surface  temperature  (SST)  cooling  (Chang  and  Anthes,  1978; 
Price,  1981;  Bender  and  Ginis,  2000;  Cione  and  Uhlhorn,  2003; 
Lin  et  al.,  2003a,b;  Tseng  et  al.,  2010).  This  SST  cooling  can  signifi¬ 
cantly  reduce  the  enthalpy  (i.e.,  sensible  heat  +  latent  heat)  flux 
supply  from  the  ocean  to  the  TC;  thus  it  is  an  important  negative 
feedback  mechanism  for  TC  intensification  (Gallacher  et  al., 
1989;  Bender  et  al.,  1993;  Emanuel,  1999;  Schade  and  Emanuel, 
1999;  Bender  and  Ginis,  2000;  Hong  et  al.,  2000;  Cione  and 
Uhlhorn,  2003;  Emanuel  et  al.,  2004;  Lin  et  al.,  2008,  2009a, 
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2013;  Balaguru  et  al.,  2012;  Moon  and  Kwon,  2012).  Over  a  thicker 
(thinner)  and  warmer  (colder)  OTS  region,  the  SST  cooling  is  gener¬ 
ally  suppressed  (enhanced)  (Price,  1981;  Shay  et  al.,  2000;  Lin 
et  al.,  2009b,  2013;  Goni  et  al.,  2009);  thus,  the  OTS  is  a  critical 
parameter  in  TC  intensity  forecasting  and  research. 

To  obtain  synoptic  OTS  observations  operationally  over  the  glo¬ 
bal  TC  basins,  in  situ  observations  are  clearly  not  the  solution.  The 
launch  of  satellite  altimeters  in  the  1990s  has  enabled  operational 
derivation  of  OTS  using  altimetry  sea  surface  height  anomaly 
(SSHA)  observations  (Goni  et  al.,  1996, 2009;  Shay  et  al.,  2000;  Wil¬ 
lis  et  al.,  2003;  Pun  et  al.,  2007,  2011;  Shay  and  Brewster,  2010). 
These  satellite-derived  OTS  fields  have  made  significant  contribu¬ 
tions  to  our  understanding  of  the  interaction  between  TCs  and 
the  ocean  (Shay  et  al.,  2000;  Goni  and  Trinanes,  2003;  Lin  et  al., 
2005,  2008;  Pun  et  al.,  2007;  Mainelli  et  al.,  2008;  Goni  et  al., 
2009;  Shay  and  Brewster,  2010).  The  existing  derivation  is  mainly 
based  on  a  two-layer  reduced  gravity  model  (hereafter  referred  to 
as  “TLM”)  where  only  two  isotherms  (i.e.,  the  depths  of  20  °C  and 


26  °C  isotherms,  or  D20  and  D26)  are  derived  to  coarsely  character¬ 
ize  the  OTS  (Shay  et  al.,  2000;  Goni  and  Trinanes,  2003;  Pun  et  al., 
2007;  Goni  et  al.,  2009;  Shay  and  Brewster,  2010).  The  TLM  has 
been  used  extensively  and  has  much  contributed  to  understanding 
the  role  ocean’s  subsurface  thermal  structure  (e.g.,  the  role  of 
warm  ocean  eddies)  plays  in  tropical  cyclone’s  intensity  change 
(Shay  et  al.,  2000;  Goni  and  Trinanes,  2003;  Mainelli  et  al.,  2008; 
Goni  et  al.,  2009;  Lin  et  al.,  2009b;  Shay  and  Brewster,  2010;  Pun 
et  al.,  2011)  and  to  monitoring  long-term  subsurface  changes  in 
typhoon  developing  region  (Pun  et  al.,  2013).  However,  improve¬ 
ment  on  vertical  resolution  of  satellite-derived  OTS  to  enhance 
ocean  characterization  is  desirable.  This  series  of  work  is  dedicated 
to  this  goal.  We  aim  to  develop  a  new  generation  of  derivation 
to  improve  ocean’s  subsurface  characterization  over  global  TC 
basins  for  better  TC  intensity  forecasting  and  research.  Since  the 
Western  North  Pacific  Ocean  (WNPO)  has  the  most  TCs  (i.e., 
typhoons)  on  the  Earth,  the  first  paper  of  this  series  (this  study) 
covers  the  WNPO. 
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Fig.  1.  (a)  Locations  of  Argo  profiles  used  in  the  present  study.  Red  dots  depict  38,556  profiles  collected  from  2000-2008,  and  blue  dots  depict  7036  profiles  from  2009-2010. 
Only  typhoon  season,  May-October,  is  considered.  Four  hydrographic  zones  are  indicated,  (b-e)  Depth-temperature  plots  of  all  the  Argo  profiles  (grey  curves)  in  September  of 
2000-2008  in  the  Gyre,  SEZ,  lEZ  and  NEZ.  In  (b-e),  the  blue  curves  depict  the  mean  profiles  and  the  blue  dots  depict  the  mean  depths  of  isotherms  from  D29  to  D4.  The 
seasonal  and  permanent  thermoclines  are  also  indicated.  (Eor  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this 
article.) 
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Table  1 

Additional  quality  control  for  Argo  profiles. 


Items 

Threshold 

Problematic  floats 

WMO^  instrument  type  852 

1  St  Points  of  temperature 

<0  °C  or  >35  °C 

Temperature  jump 

>1  °C 

Maximum  depth 

<200  m 

Depth  jump 

>300  m 

Valid  points 

<10  points 

^  WMO:  World  Meteorological  Organization. 


Fig.  la  depicts  the  WNPO  domain  (120-170°E  and  0-40°N). 
According  to  Lin  et  al.  (2005,  2008)  and  Pun  et  al.  (2007),  four 
zones  with  distinctive  TC-related  OTS  characteristics  can  be  found 
in  the  WNPO.  These  zones  are  the  Gyre,  Southern  Eddy-rich  Zone 
(SEZ),  Intermediate  Eddy  Zone  (lEZ),  and  Northern  Eddy-rich  Zone 
(NEZ).  The  Gyre  (125-170°E  and  10-21°N)  is  the  region  to  the 
south  where  the  deepest  and  warmest  OTS  is  located.  In  the  SEZ 
(127-170°E  and  21-26°N)  and  NEZ  (140-170°E  and  30-40°N), 
warm  and  cold  ocean  eddies  frequently  appear.  Thus,  in  the  SEZ 
and  NEZ,  the  OTS  can  greatly  fluctuate.  Lastly,  the  lEZ  (135- 
170°E  and  26-30°N)  is  the  intermediate  region,  in  between  the 
two  eddy-rich  zones  of  SEZ  and  NEZ  (Pun  et  al.,  2007).  Prom  the 
perspective  of  TC,  the  OTS  shows  distinctive  differences  among 
these  zones  (Fig.  Ib-e).  The  Gyre  region  has  the  warmest  subsur¬ 
face  temperature,  as  characterized  by  the  deepest  D20  and  D26. 
Moving  north  to  the  SEZ,  lEZ,  and  NEZ,  the  OTS  shoals  significantly 
poleward.  In  the  NEZ,  the  OTS  is  the  shallowest  and  the  mean  D20 
is  only  around  60  m,  much  shallower  compared  to  the  D20  of 
~200m  in  the  Gyre  region.  Purther  details  of  these  four  zones 
and  the  hydrographical  setting  of  the  WNPO  can  be  found  in  the 
Supplementary  Online  Material  (SOM). 

Pun  et  al.  (2007)  assessed  the  accuracy  of  the  existing  TLM  der¬ 
ivation  in  the  WNPO.  It  was  found  that  D20  can  be  grossly  overes¬ 
timated  by  more  than  100%  in  the  northern  WNPO  (i.e.,  the  NEZ 
region).  Therefore,  in  addition  to  increase  the  vertical  resolution, 
we  intend  to  find  a  new  method  to  improve  the  OTS  derivation 
in  the  WNPO. 

Derivation  based  on  linear  regression  is  an  approach  which  has 
been  used  in  climate-related  studies  (Cheney,  1982;  Carnes  et  al., 
1990;  White  and  Tai,  1995;  Gilson  et  al.,  1998;  Willis  et  al.,  2003, 
2004)  but  not  in  TC-related  studies.  This  is  because  meaningful  lin¬ 
ear  regressions  can  only  be  established  with  a  large  number  of 
in  situ  profiles  and  there  was  not  enough  data  for  TC-related 


derivation.  To  date,  with  the  great  enhancement  in  in  situ  observa¬ 
tions  from  the  Argo  floats  (Gould  et  al.,  2004;  Roemmich  et  al., 
2004),  it  has  become  possible  to  apply  regression  method  to  derive 
the  OTS  for  TC-related  purposes  (Guinehut  et  al.,  2012).  In  this 
study,  we  derive  daily,  0.25°  OTS  from  satellite  altimetry  over  the 
WNPO  during  the  typhoon  season  (May-October)  using  a  regres¬ 
sion  method.  This  is  established  based  on  all  available  (>38,000) 
in  situ  Argo  float  profiles  collected  during  the  typhoon  season  in 
2000-2008. 

2.  Data 

2.1.  Argo  profiles 

Argo  is  an  autonomous  float  that  continuously  measures  tem¬ 
perature  and  salinity  profiles  in  the  ocean  (typically  from  surface 
to  2000  m  depth)  at  a  10-day  interval  (Gould  et  al.,  2004;  Roem¬ 
mich  et  al.,  2004).  Currently,  over  3500  Argo  floats  are  collecting 
data  at  any  given  time  period,  providing  an  unprecedented  number 
of  real-time  and  high-quality  in  situ  measurements  over  the  global 
oceans.  In  this  study,  1 1  years  of  Argo  thermal  profiles  from  2000- 
2010  are  used  (Pig.  la).  Profiles  in  the  first  nine  years  of  2000-2008 
are  used  to  build  the  linear  regressions  between  satellite  SSHA  and 
subsurface  temperature  variation,  whereas  the  profiles  in  the  last 
two  years  of  2009  and  2010  are  used  for  the  validation  of  the  de¬ 
rived  OTS. 

Although  the  profiles  provided  by  the  Argo  data  center 
(www.argo.net)  have  already  gone  through  a  standard  quality  con¬ 
trol  (QC)  procedure,  additional  QC  is  performed  to  ensure  that  all 
the  profiles  used  in  this  study  are  in  high  quality.  According  to  Wil¬ 
lis  et  al.  (2009),  cold  biases  were  found  in  some  types  of  Argo  floats 
deployed  between  mid-2003  and  2006.  for  this  reason,  the  profiles 
from  the  problematic  floats  are  excluded  in  this  study.  All  of  the 
items  and  their  thresholds  for  the  additional  QC  are  listed  in  Ta¬ 
ble  1.  On  average,  about  10%  of  the  profiles  were  ruled  out  by  this 
additional  QC. 

2.2.  Satellite  SSHA  and  SST 

The  near  real-time  (NRT)  daily  and  delay-time  (DT)  7-day 
merged  and  gridded  SSHA  products  were  employed  in  this  study. 
Both  products  have  the  spatial  resolution  of  0.25°.  They  are  opera¬ 
tionally  generated  by  the  Segment  Sol  multimissions  d’ALTimtoie, 
d’Orbitographie  et  de  localisation  precise/Data  Unification  and 
Altimeter  Combination  System  (SSALTO/DUACS)  and  distributed 


Fig.  2.  The  monthly  correlation  (R)  between  SSHA  and  the  displacements  of  the  isotherms  from  29  °C  to  4  °C,  in  (a)  the  Gyre,  (b)  the  SEZ,  (c)  the  lEZ  and  (d)  the  NEZ.  The  color 
code  for  each  curve  is  shown  in  (a).  The  dashed  lines  represent  the  99%  significance  level,  and  the  dotted  lines  depict  R  =  0.6.  (For  interpretation  of  the  references  to  color  in 
this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  3.  Schematic  diagram  shows  how  a  D20  regression  is  developed,  (a)  An  exemplary  depth-temperature  profile  pair  of  an  Argo  in  situ  profile  (black)  and  WOAOl 
climatological  profile  (blue).  The  left  axis  is  depth  and  the  right  axis  is  SSHA.  The  D20  displacement,  i.e.  AD20,  is  the  difference  between  the  in  situ  and  climatology,  and  its 
corresponding  SSHA  value  is  indicated  as  a  yellow  vertical  bar.  (b)  The  scatter  plot  of  AD20  and  SSHA  consists  of  a  number  of  observation  pairs.  Based  on  these  pairs,  a  linear 
regression  is  developed. 


via  the  Archiving,  Validation  and  Interpretation  of  Satellite 
Oceanographic  data  (AVISO)  web  site  (www.aviso.oceanobs.com). 
In  these  products,  SSHA  measurements  from  all  available  altimetry 
missions  are  merged  to  maximize  the  observations  (Ducet  et  al., 
2000;  Pascual  et  al.,  2006,  2009).  These  missions  include  Ocean 
TOPography  Experiment  (TOPEX)/Poseidon  (T/P),  Jason-1,  Jason-2, 
Geosat  Eollow-On  (GEO),  Envisat,  and  European  Remote  sensing 
Satellite  (ERS)  series  (Eu  et  al.,  2003;  Le  Traon  and  Dibarboure, 
2004).  This  AVISO  merged  dataset  is  the  most  well-known  and 
was  widely  used  in  many  studies  (Niiler  et  al.,  2003;  Willis  et  al., 
2004;  Qiu  and  Chen,  2010;  Cazenave  and  Remy,  2011;  Chelton 
et  al.,  201 1 ).  To  match  the  Argo  observations,  the  SSHA  data  be¬ 
tween  2000  and  2010  is  used.  Initially,  the  7-day  DT  product  is 
used.  When  the  higher  temporal-resolution  (daily)  data  from  the 
NRT  dataset  became  available  after  5  July  2006,  the  daily,  NRT  data 
is  used  from  5  July  2006  to  2010. 

Besides  the  SSHA,  satellite-observed  SST  is  also  used.  The  ver¬ 
sion-3  fusion  SST  is  used  in  this  study,  which  is  from  the  cloud  pe¬ 
netrating  Tropical  Rainfall  Measuring  Mission/Microwave  Imager 
(TRMM/TMI)  and  Advanced  Microwave  Sounding  Radiometer  for 
the  Earth  Observing  System  (AMSR-E).  This  dataset  is  operationally 
provided  by  the  Remote  Sensing  Systems  with  a  nominal  accuracy 
of  ~0.5  °C  (Wentz  et  al.,  2000).  The  spatial  and  temporal  resolu¬ 
tions  are  0.25°  and  daily,  respectively. 


3.  Relationship  between  SSHA  and  depth  variation  of  isotherms 

Based  on  the  assumption  that  the  change  in  sea  surface  height  is 
primarily  contributed  by  the  steric  effect  (i.e.,  thermal  expansion 
and  saline  contraction;  Gill  and  Niiler,  1973),  the  variations  in  sub¬ 
surface  properties,  such  as  temperature  and  salinity,  are  closely  re¬ 
lated  to  the  variation  in  SSHA,  which  is  now  accurately  detected  by 
space-based  altimeters.  In  the  tropical  and  subtropical  oceans,  the 
changes  in  sea  surface  height  are  largely  dominated  by  tempera¬ 
ture  changes  in  the  water  column  due  to  thermal  expansion  (White 
and  Tai,  1995;  Gilson  et  al.,  1998;  Willis  et  al.,  2003).  Thus,  positive 
SSHAs  generally  reflect  subsurface  temperatures  that  are  warmer 


than  the  climatology  (isotherm  deepening  below  the  mean  level; 
Eig.  3a).  In  contrast,  negative  SSHAs  represent  cooler  temperatures. 

To  start,  we  calculate  monthly  correlation  coefficients  between 
SSHA  and  the  depth  variation  in  each  isotherm  from  D4  to  D29 
(Eig.  2).  The  variation  of  the  isotherm  depth  is  the  vertical  displace¬ 
ment  of  the  in  situ  isotherm  (i.e.,  from  Argo)  with  respect  to  the 
monthly  climatology  held  of  the  World  Ocean  Atlas  2001 
(WOAOl;  Stephens  et  al.,  2002).  In  other  words,  the  displacement 
is  the  difference  between  Argo  and  WOAOl  climatology  within 
the  same  month  (Eig.  3a).  Thus  the  correlations  in  Eig.  2  are  based 
on  the  observations  (i.e.,  satellite  SSHA  and  Argo  temperature  pro¬ 
files  from  2000  to  2008).  It  can  be  seen  that  the  correlation  varies 
from  isotherm  to  isotherm,  from  region  to  region,  and  from  month 
to  month.  In  the  Gyre  (Eig.  2a),  strong  correlation  of  >0.6  can  be 
found  in  most  isotherms  between  D12  to  D27  with  a  maximum 
(~0.75)  found  between  D18  and  D22.  The  correlation  also  varies 
with  month,  with  slightly  lower  correlations  in  May/June  and  a 
higher  correlation  in  October.  In  the  SEZ,  high  correlation  of  >0.6 
is  found  almost  in  the  entire  column  of  the  isotherms  (i.e.,  D4- 
D27),  with  a  maximum  (~0.8)  found  between  D19  and  D23.  Inter¬ 
estingly,  a  local  minimum  is  found  in  D17.  This  may  be  due  to  the 
fact  that  D1 7  is  located  in  the  transition  layer  between  the  seasonal 
and  permanent  thermoclines  and  is  expected  to  be  less  correlated 
to  SSHA  (Eig.  Ic). 

Moving  north,  the  situation  in  the  lEZ  becomes  more  complex; 
the  isotherms  with  a  higher  correlation  of  >0.6  are  generally 
located  below  D17  (Eig.  2c),  all  within  the  permanent  thermocline 
(Eig.  Id).  In  contrast,  the  isotherms  above  D17  belong  to  the 
seasonal  thermocline  and  have  lower  correlations  (temporal  max¬ 
imum  correlations  are  found  in  August/September).  Similar  to  the 
SEZ,  a  local  minimum  in  the  lEZ  is  also  found  in  D17,  indicating 
that  the  transition  layer  is  also  around  this  isotherm.  Einally,  in 
the  NEZ  (Eig.  2d),  the  isotherms  with  the  highest  correlation 
(~0.8-0.9)  are  found  below  D17  (i.e.,  the  layer  of  the  permanent 
thermocline;  Eig.  le). 

Erom  above,  it  can  be  seen  that  though  SSHA  is  only  used  to  re¬ 
trieve  D20  and  D26  in  the  current  TLM  (Shay  et  al.,  2000;  Goni  and 
Trinanes,  2003;  Pun  et  al.,  2007;  Goni  et  al.,  2009;  Shay  and 
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Fig.  4.  The  spatial  distribution  of  (a)  the  slope  and  (b)  the  offset  of  the  Eq.  (2)  for  D20  in  September,  (c)  The  spatial  distribution  of  the  correlation  between  D20  and  SSHA.  (d) 
The  number  of  profiles  used  for  developing  the  regression  at  each  0.25°  x  0.25°  grid  point.  Grey  areas  represent  the  significance  level  <99%. 


Brewster,  2010),  other  isotherms  in  the  WNPO  can  have  similar  or 
even  stronger  correlations.  This  suggests  that  it  is  possible  to  re¬ 
trieve  other  isotherms  using  SSHA  and  not  limited  to  only  D20 
and  D26. 

4.  Linear  regression 

A  total  of  38,556  Argo  profiles  are  used  to  develop  the  regres¬ 
sions  between  SSHA  and  subsurface  isotherm  depth  variation.  We 
construct  a  regression  at  every  0.25°  x  0.25°  grid  and  on  a 
monthly  basis.  The  vertical  resolution  is  at  1  °C  interval  from 
D4  to  D29  (26  layers).  To  have  a  robust  and  representative  regres¬ 
sion,  a  large  number  of  in  situ  profiles  are  required.  However,  in  a 
fine  resolution  grid  of  0.25°,  it  is  unlikely  to  have  a  significant 
number  of  in  situ  profiles.  To  solve  this  problem,  we  use  a  mov¬ 
ing-window  approach.  For  each  grid,  profiles  within  a  window 
of  10°  longitude  by  5°  latitude  centered  at  the  grid  are  binned 
to  calculate  the  regressions.  The  reason  for  using  10°  longitude 
is  because  there  is  less  variability  in  the  east-west  direction 
(Fig.  1).  The  regressions  are  constructed  by  linearly  regressing 
the  isotherm  displacements  onto  their  co-incident  and  co-located 
SSHA  values.  The  displacement  thus  is  the  isotherm  depth  differ¬ 
ence  between  in  situ  and  climatological  profiles,  representing  the 
deviation  from  the  climatology  (Fig.  3a).  Here,  the  climatology 
field  is  from  the  monthly,  0.25°  temperature  data  of  the  WOAOl 
(Stephens  et  al.,  2002).  The  isotherm  displacement  (AD)  in  each 
in  situ  profile  is  derived  as: 

ADi  =  Di-N,  1  =  29, 28,..., 4,  (1) 


where  D  and  D  are  the  in  situ  and  climatological  depths  of  isotherm, 
respectively;  and  i  is  the  isotherm  from  29  °C  to  4  °C. 

The  displacements  in  each  isotherm  of  the  in  situ  profiles  are 
linearly  regressed  onto  their  corresponding  SSHA  by  least  square 
fit  (Fig.  3b).  The  monthly  regression  equation  for  each  isotherm 
at  each  grid  is  obtained  using: 

ADi  =  Cl(x,y,m,i)  ^  tj  b(^x.y,m,i)^  (2) 

where  r]  is  the  satellite-observed  SSHA;  a  and  b  are  the  linear 
regression  coefficients  representing  the  slope  and  offset,  respec¬ 
tively;  both  coefficients  are  a  function  of  location  {x,y\  interval  of 
0.25°),  isotherm  (i;  D29-D4),  and  month  (m;  May-October).  With 
this  set  of  regression  equations,  the  displacement  of  an  isotherm 
(ADi)  can  be  estimated  from  the  SSHA.  The  actual  depth  of  the 
isotherm  (Di)  is  then  retrieved  by  adding  the  corresponding  clima¬ 
tological  depth.  To  precisely  demonstrate  the  derivation  of  the 
regression  equations  (Eq.  (2))  from  Argo  profiles  and  satellite  SSHA, 
we  use  D20  as  an  example  to  make  a  schematic  diagram  (Fig.  3).  An 
Argo  profile  in  Fig.  3a  represents  one  arbitrary  point  in  Fig.  3b. 
Firstly,  AD20  is  obtained  by  subtracting  the  climatological  D20  from 
in  situ  D20,  and  then  paired  with  co-located  and  co-incident  SSHA 
(yellow  bar  in  Fig.  3a).  Given  a  large  number  of  AD20-SSHA  pairs 
(red  dots  in  Fig.  3b),  a  linear  regression  for  D20  is  derived.  Equations 
for  other  isotherms  are  derived  in  the  same  way. 

Eig.  4  presents  the  regression  coefficients  of  Eq.  (2)  for  D20  in 
September.  The  regression  is  developed  based  on  ~50-300  profiles 
for  each  grid  (Eig.  4d).  One  can  see  that  the  slope  (a)  of  the  regres¬ 
sions  exhibits  spatial  variability.  The  steepest  slope  is  found  near 
15-25°N,  about  2.5  m  cm“\  meaning  that  1cm  of  SSHA  can 
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Fig.  5.  (a)  The  isotherm-SSHA  correlation  from  the  equator  to  40°N  along  145°E  in  September,  (b)  Same  as  in  (a),  but  for  the  slope  of  the  regression.  The  slope  less  than  99% 
significance  is  marked  in  grey.  The  latitudinal  band  for  the  Gyre,  SEZ,  lEZ  and  NEZ  are  indicated. 


produce  2.5  m  variation  in  D20.  The  offset  {b)  is  generally  negative 
(Fig.  4b).  Fig.  4c  shows  the  correlation  in  September.  Basically,  the 
correlation  between  SSHA  and  D20  is  high  (~0.7).  However,  low 
correlation  can  be  found  in  local  regions  to  the  north. 

Fig.  5a  illustrates  the  isotherm-SSHA  correlation  from  the 
equator  to  40°N  along  145°E  in  September,  to  demonstrate  the 
relationship  in  the  vertical  direction.  It  can  be  found  that  the  cor¬ 
relation  pattern  varies  significantly  with  isotherm  and  latitude. 
Basically,  the  correlation  pattern  is  roughly  consistent  with  the 
four  zones  (Fig.  1),  but  small-scale  spatial  variability  still  exists. 
Fig.  5b  shows  the  corresponding  cross  section  for  the  slope  of 
the  regression.  The  steepest  slope  is  found  in  the  lower  isotherms 
(i.e.,  D4-D17)  in  the  lEZ  and  NEZ,  in  which  the  slope  is  up  to  ~4- 
5  m  cm“\  It  also  corresponds  to  the  high  correlation  area  found 
in  Eig.  5a.  In  the  Gyre  and  SEZ,  the  larger  slope  of  ~2.5  m  cm“^ 
is  found  in  the  upper  isotherms  (typically  ~D25-D19).  The  grey 
regions  in  Eig.  5b  are  locations  where  slope  does  not  meet  the 
99%  statistical  significance  level  and  no  retrieval  on  these  regions 
is  performed. 

After  obtaining  the  26  isotherms,  the  OTS  derivation  is  com¬ 
pleted  through  inclusion  of  the  satellite-observed  SST  (Section  2) 
and  climatological  monthly  mixed  layer  depth  (MED)  from  the 
U.S.  Naval  Research  Laboratory  (NRL;  Kara  et  al.,  2002).  Thus,  the 
OTS  from  surface  to  the  depth  of  ^lOOOrn  (i.e.,  D4)  can  be  ob¬ 
tained,  and  can  be  compared  to  the  two  isotherms  based  on  the 


existing  method  (Eig.  6).  Hereafter,  this  regression  method  is  re¬ 
ferred  to  as  the  “REGWNP”  (regression  for  the  western  North  Paci¬ 
fic).  It  should  be  noted  that  the  derived  OTS  is  regular  temperature, 
not  potential  temperature.  Given  the  daily  SSHA  and  SST  maps,  the 
corresponding  daily  OTS  can  be  derived  through  the  monthly 
regression,  climatological  isotherm  depths  and  MLD  within  the 
same  month.  Pig.  7  displays  exemplary  profile  comparison  in  Sep¬ 
tember  2010  from  the  REGWNP  (green  curves),  Argo  (black  curves) 
and  TLM  (red  curves)  in  the  four  zones.  It  can  be  seen  that  the  REG¬ 
WNP  profiles  are  generally  much  closer  to  the  in  situ  Argo  profiles, 
as  compared  to  those  from  the  TLM.  Detailed  validation  results  are 
presented  in  Section  5. 

5.  Validation  and  comparison  with  the  TLM 

An  independent  set  of  7036  Argo  profiles  from  the  2009  and 
2010  typhoon  seasons  were  used  to  validate  the  regression-de¬ 
rived  OTS.  The  locations  of  these  profiles  are  shown  in  Pig.  la  as 
blue  dots.  Series  of  validation  for  the  REGWNP  were  conducted, 
and  comparison  with  the  existing  TLM  was  made.  Root-mean- 
square  differences  (RMSDs),  biases,  and  errors  were  calculated  be¬ 
tween  the  observed  (i.e.,  Argo  profiles)  and  the  REGWNP-derived 
isotherms  to  quantify  the  performance  of  the  REGWNP.  Here,  error 
is  defined  as  the  percentage  of  RMSD  in  the  observed  mean: 


Fig.  6.  Schematic  of  (a)  REGWNP-derived  profile  and  (b)  TLM-derived  profile,  (c)  Zoom-in  area  of  (a)  compares  the  profiles  from  the  REGWNP  (green)  and  Argo  (black),  (d) 
Zoom-in  area  of  (b)  compares  the  profiles  from  the  TLM  (red)  and  Argo  (black).  Shaded  areas  in  (c)  and  (d)  indicate  the  UOHC,  their  values  are  also  shown.  (For  interpretation 
of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  7.  Comparison  between  the  profiles  from  the  REGWNP  (green),  TLM  (red)  and  Argo  (black).  From  the  upper  to  lower  rows  represent  the  profiles  in  the  NEZ,  lEZ,  SEZ  and 
Gyre.  Each  zone  has  four  randomly-selected  profiles  in  September  of  2010.  The  TLM  profile  is  not  available  in  the  upper  right  plot.  The  horizontal  and  vertical  dashed  lines  depict 
26  °C  and  the  depth  of  100  m,  respectively.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  8.  Scatter  plots  of  the  REGWNP  (green)  and  TLM  (red)  against  Argo  observations  for  (a)  D8,  (b)  D14,  (c)  D20,  (d)  D26,  (e)  UOHC,  and  (f)  TlOO  for  all  four  zones  (the  Gyre, 
SEZ,  lEZ  and  NEZ).  The  Argo  profiles  compared  are  from  May-October  of  2009-2010  (blue  dots  in  Fig.  la).  Sample  number,  RMSD,  correlation  and  bias  with  respect  to  Argo  are 
shown.  The  TLM  cannot  derive  D8  in  (a)  and  D14  in  (b).  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this 
article.) 
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In  addition  to  the  derived  isotherms,  two  other  typhoon-inten¬ 
sification-related  parameters  calculated  from  the  isotherms  were 
also  calculated  and  validated.  These  parameters  are  TlOO  (top 
100-m  averaged  ocean  temperature;  Price,  2009)  and  tropical  cy¬ 
clone  heat  potential  [TCHP;  also  called  upper  ocean  heat  content 
(UOHC)]  (Leipper  and  Volgenau,  1972;  Shay  et  al.,  2000;  Goni 
and  Trinanes,  2003;  Goni  et  al.,  2009;  Shay  and  Brewster,  2010). 
From  Shay  et  al.  (2000)  and  Goni  and  Trinanes  (2003),  TCHP  (or 
UOHC)  is  the  heat  content  exceeding  26  °C  in  the  ocean.  It  can  be 
calculated  from  the  derived-OTS  profile  using: 


where  Cp  is  the  heat  capacity  of  the  water  at  the  constant  pressure  of 
4178  J  kg“^°C“\  p  is  the  average  density  of  the  upper  ocean  taken 
as  1026  kg  m“^,  and  AT  is  the  difference  between  r(z)  and  26  °C 
at  the  depth  z. 

For  the  comparison  with  the  TLM,  OTS  is  also  calculated  by  the 
TLM.  First,  D20  is  derived  from  SSHA  thought  the  two-layer  re¬ 
duced  gravity  scheme  proposed  by  Goni  et  al.  (1996)  and  Shay 
et  al.  (2000): 

D20  =  D20  +  w,  (5) 
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Fig.  9.  (a)  Overall  RMSD  of  REGWNP-derived  isotherms  from  D29-D4  over  all  four  zones  (the  Gyre,  SEZ,  lEZ  and  NEZ)  during  May-October  of  2009-2010.  (b)  Same  as  (a)  but 
for  observed  mean  depth  of  isotherms  from  Argo,  (c)  Same  as  (a),  but  for  error  of  REGWNP-derived  isotherms,  (d)  Same  as  (a),  but  for  correlation  between  REGWNP-derived 
and  Argo  isotherms. 
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Fig.  10.  Regional  D20  scatter  plots  of  the  REGWNP  (green)  and  TLM  (red)  against  Argo  observations  for  the  (a)  Gyre,  (b)  SEZ,  (c)  lEZ,  and  (d)  NEZ.  The  Argo  profiles  compared 
are  from  May-October  of  2009-2010  (blue  dots  in  Fig.  la).  Sample  number,  RMSD,  correlation  and  bias  with  respect  to  Argo  are  shown.  (For  interpretation  of  the  references  to 
color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


where  D20,  pi  and  p2  are  the  climatological  D20,  upper  layer  and 
lower  layer  densities,  respectively,  which  are  obtained  from  the 
WOAOl  monthly  climatology  fields.  After  obtaining  D20,  D26  is  de¬ 
rived  by  multiplying  a  climatological  ratio  (D26/D20)  from  the 
WOAOl.  Finally,  OTS  is  obtained  with  the  satellite  SST  and  NRL’s 
monthly  climatological  MLD  (Fig.  6b).  More  detailed  information 
about  TLM  derivation  in  the  WNPO  can  be  found  in  Pun  et  al. 
(2007). 

5.1.  General  validation 

Fig.  8  shows  that  the  REGWNP-derived  isotherms  (green  dots) 
are  reasonably  constrained  in  RMSD  with  no  evident  biases.  For 
the  two  well-known  typhoon-related  isotherms  used  in  the  TLM 
(Leipper  and  Volgenau,  1972;  Shay  et  al.,  2000;  Goni  and  Trin- 
anes,  2003;  Lin  et  al.,  2008),  i.e.,  D26  and  D20,  the  RMSDs  (their 
biases)  are  15  m  (-2  m)  and  20  m  (-3  m),  respectively,  as  shown 
in  Fig.  8c  and  d.  In  comparison,  there  are  significantly  larger 


biases  (i.e.,  over-estimation)  derived  from  the  existing  TLM  for 
D20  (red  dots  in  Fig.  8c),  which  was  also  pointed  out  by  Pun 
et  al.  (2007).  For  the  D20  from  the  TLM,  the  RMSD  is  as  larger 
as  42  m  (i.e.,  twice  the  RMSD  of  the  REGWNP).  The  performance 
of  D26  from  the  TLM  is  better  than  that  of  D20,  though  its  RMSD 
is  Still  larger  than  those  from  the  REGWNP  (Pig.  8d).  In  terms  of 
TCHP  (or  UOHC)  and  TlOO,  it  can  be  seen  that  though  the  REG¬ 
WNP  performs  better  than  the  TLM,  both  can  obtain  reasonable 
derivation,  with  RMSD  being  ~15-16kjcm“^  for  UOHC  and 
-0.6-0.8  °C  for  TlOO  (Pig.  8e  and  f). 

Because  the  REGWNP  can  derive  26  isotherms,  we  validate  the 
rest  of  the  isotherms  as  well.  Pig.  8a  and  b  depict  the  scatter  plots 
for  the  two  deeper  isotherms,  i.e.,  D8  and  D14,  while  the  complete 
validation  for  all  26  isotherms  are  summarized  in  Pig.  9.  As  in 
Pig.  8a  and  b,  the  RMSD  for  the  lower  isotherms  are  also  reasonably 
constrained  with  most  dots  fall  near  the  1:1  diagonal  line  between 
the  observed  and  the  derived  values.  Por  D14  (D8),  the  RMSD  is 
29  m  (38  m).  Because  these  are  deeper  isotherms,  located  at 
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Fig.  11.  Monthly  RMSD  (1st  row),  bias  (2nd  row),  error  percentage  (3rd  row),  and  sample  number  (4th  row)  for  the  REGWNP-derived  D20  (1st  column  from  the  left),  D26 
(2nd  column),  UOHC  (3rd  column),  and  TlOO  (4th  column).  Green,  blue,  orange  and  red  curves  represent  the  Gyre,  SEZ,  lEZ  and  NEZ,  respectively.  Dashed  and  solid  curves 
represent  the  TLM  and  REGWNP,  respectively.  The  units  for  RMSD,  bias,  and  error  are  shown.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  this  article.) 


~400  m  and  900  m  (Fig.  9b),  the  errors  (in  terms  of  RMSD/observed 
mean)  are  actually  smaller  (Fig.  9c).  It  can  be  seen  that  for  the 
deeper  isotherms  of  D18  and  below,  the  errors  can  be  constrained 
to  within  5-10%.  For  the  upper  ocean,  the  errors  are  larger,  though 
they  are  still  within  20%  (Fig.  9c). 

The  quality  of  the  REGWNP-derived  isotherms  can  also  be  seen 
from  Fig.  9d,  which  depicts  the  correlation  between  the  derived 
and  the  observed  isotherms.  High  correlation  of  ~0.9  exists  for 
all  isotherms  but  the  very  top  isotherms  of  D28  and  D29  near 
the  mixed  layer.  Fig.  9d  also  illustrates  one  important  point,  that 
is,  though  traditionally  only  D20  and  D26  are  derived  using  satel¬ 
lite  altimetry  (Shay  et  al.,  2000;  Goni  and  Trinanes,  2003),  many 
other  isotherms  can  also  be  derived  with  as  good,  or  even  better, 
accuracy.  This  is  encouraging  because  high-resolution  and  much 
more  complete  subsurface  information  can  be  derived  from  the 
REGWNP.  As  will  be  shown  in  Section  6,  this  improvement  can 
lead  to  more  accurate  estimation  of  air-sea  enthalpy  (latent 
heat  +  sensible  heat)  flux  and  typhoon-induced  ocean  cooling, 
both  are  critical  quantities  for  research  that  concerns  typhoon 
intensiflcation  and  typhoon-ocean  interaction  (Price,  1981, 
2009;  Emanuel,  1999;  Bender  and  Ginis,  2000;  Clone  and  Uhl- 
horn,  2003;  Emanuel  et  al.,  2004;  Lin  et  al.,  2005,  2009a,  2013; 
Clone  et  al.,  2013). 

Eurther,  it  can  be  seen  that  the  isotherm  derivation  is  not  re¬ 
stricted  to  the  upper  ocean,  but  can  be  extended  to  deeper  ocean 
(i.e.,  D4,  typically  ~1000  m;  see  Eigs.  Ib-le  and  9b).  This  suggests 
that  the  ocean  thermal  profiles  derived  from  the  REGWNP  can 
also  be  applied  to  other  scientific  issues  in  oceanography  over 
the  WNPO,  not  just  in  the  held  of  typhoon-ocean  interaction  that 
typically  takes  place  in  the  top  300  m  (Price,  1981,  2009;  Lin 
et  al.,  2005;  D’Asaro  et  al.,  2011;  Lin,  2012;  Mrvaljevic  et  aL, 
2013). 


5.2.  Regional  and  monthly  validations 

In  addition  to  general  validation,  regional  validation  is  per¬ 
formed  to  assess  the  accuracy  in  different  parts  of  the  WNPO 
(Pig.  la).  This  provides  additional  information  on  the  regional 
dependence  of  the  performance;  especially,  it  was  noted  by  Pun 
et  al.  (2007)  that  the  D20  from  the  TLM  suffers  large  over-estima¬ 
tion  over  the  NEZ  region.  Pig.  10  depicts  the  results  of  D20  deriva¬ 
tion.  It  can  be  seen  that  for  the  two  zones  at  the  southern  part  of 
the  WNPO,  i.e..  Gyre  and  SEZ,  the  REGWNP  is  well  constrained  with 
RMSD  of  16  m  and  20  m,  respectively  (green  dots  in  Pig.  10a  and  b). 
Towards  the  northern  part  of  the  WNPO,  i.e.,  the  lEZ  and  the  NEZ 
regions,  the  RMSD  increases  to  22  m  and  25  m  and  the  correlation 
of  ~0.7  is  also  lower  than  that  in  the  Gyre  and  the  SEZ  regions 
(with  correlation  ~0.83-0.84).  This  shows  that  the  performance 
of  the  REGWNP  is  better  over  the  southern  WNPO  (south  of 
26°N)  than  over  the  northern  part  (Pig.  la).  As  reported  by  Pun 
et  al.  (2007)  and  in  the  SOM,  the  northern  regions  are  more  com¬ 
plex  oceanographically,  including  regions  where  the  Oyashio  and 
the  Kuroshio  meet  and  the  Kuroshio  meanders  (Yasuda  et  al., 
1992;  Qiu,  1999).  It  is  therefore  understandable  that  satellite  SSHA 
derivation  does  not  perform  as  well  as  it  does  in  the  southern 
regions. 

The  same  situation  is  found  in  the  TLM,  but  the  northern  degra¬ 
dation  in  performance  is  much  more  severe  (red  dots  in  Pig.  10c 
and  d).  It  can  be  seen  that  the  RMSD  of  83  m  over  the  NEZ  is  more 
than  300%  as  compared  to  the  REGWNP  method  (RMSD  25  m) 
(Pig.  lOd).  Also,  very  large  biases  (over-estimation)  of  46  m  are 
found  in  the  TLM,  as  compare  to  the  small  bias  of  -3  m  in  the  REG¬ 
WNP.  Using  the  TLM,  the  correlation  between  the  observed  and  the 
derived  D20  drops  to  0.47,  as  compared  to  the  higher  correlation  of 
0.71  found  using  the  REGWNP  (Fig.  lOd).  As  quantifled  in  Fig.  11c, 
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Fig.  12.  Same  as  in  Fig.  10,  but  for  TIOO. 


this  is  equivalent  to  more  than  100%  error  (red  dashed  curve),  as 
compared  to  the  ~30%  error  (red  solid  curve)  for  the  REGWNP 
method.  Over  the  two  southern  regions,  the  performance  of  the 
TLM  is  much  better,  though  it  is  still  not  as  good  as  that  of  the  REG¬ 
WNP  (Fig.  10a  and  b). 

The  above  comparison  shows  the  regional  dependence  of  the 
performance  in  D20  derivation.  Over  the  southern  regions,  both 
the  REGWNP  and  the  TLM  are  usable,  though  the  REGWNP  has 
higher  accuracy.  Over  the  northern  region,  the  REGWNP  has  less 
accuracy  and  should  be  used  with  caution  (Figs.  10  and  11).  It  is 
not  advisable  to  use  the  TLM  over  the  two  northern  regions,  given 
the  low  performance.  Similar  situation  is  also  found  in  TIOO  esti¬ 
mation  (Fig.  12). 

Quantification  of  the  regional  performance  for  D20,  D26,  UOHC, 
and  TIOO  using  both  methods  is  given  in  Fig.  11.  Also,  temporal 
performance  is  assessed.  In  general,  there  is  less  accuracy  in  the 
earlier  typhoon  season  (i.e..  May  and  June)  for  both  methods. 


And  the  best  performance  for  both  methods  is  found  in  July-Sep- 
tember  (3rd  row  in  Fig.  11). 


6.  Discussion 

6.1.  Applicability  of  REGWNP- derived  typhoon-related  parameters 

According  to  the  validation  (Section  5),  it  is  evident  that  the 
REGWNP  can  reasonably  estimate  the  typhoon-related  parame¬ 
ters:  D20,  D26,  UOHC  and  TIOO  in  the  WNPO.  As  a  convenient 
guide,  we  experimentally  took  35%  of  error  as  a  threshold  to  deter¬ 
mine  whether  the  parameters  derived  by  the  REGWNP  are  of 
acceptable  accuracy.  Fig.  13  shows  the  applicable  zones  and 
months  of  the  REGWNP  for  each  typhoon-related  parameter.  It 
can  be  seen  that  REGWNP  is  well  applicable  in  the  Gyre  throughout 
the  May-October  typhoon  season.  Similar  applicability  is  also 
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Fig.  13.  Applicability  chart  for  REGWNP-derived  D20,  D26,  UOHC  and  TlOO  for  typhoon  applications.  From  upper  to  lower  rows  are  for  the  NEZ,  lEZ,  SEZ  and  Gyre, 
respectively.  The  months,  May-October,  are  indicated  in  each  pie  chart.  The  applicability  threshold  is  defined  as  an  error  of  35%. 


found  in  the  SEZ,  except  some  part  of  the  early  season  (i.e.,  May 
and  June).  Due  to  colder  OTS  in  the  earlier  typhoon  season  (May 
and  June)  in  the  two  northern  zones  (i.e.,  NEZ  and  lEZ),  the  error 
is  larger  because  the  mean  values  are  usually  small  (see  Eq.  (3)). 
Therefore,  for  these  northern  zones,  it  is  more  difficult  to  achieve 
<35%  error. 

6.2.  Improvement  on  SST  cooling  and  enthalpy  flux  estimations 

When  considering  ocean’s  control  over  typhoon  intensity 
change,  the  most  important  quantities  to  be  obtained  are  the  dur- 
ing-typhoon  SST  cooling  and  air-sea  enthalpy  (latent  heat  +  sensi¬ 
ble  heat)  fluxes  (Emanuel,  1999;  Bender  and  Ginis,  2000;  Gone  and 
Uhlhorn,  2003;  Emanuel  et  al.,  2004;  Lin  et  al.,  2005,  2009a,  2013; 
Gone  et  al.,  2013).  Though  in  situ  sampling  has  been  much  im¬ 
proved  through  the  Argo  float  deployments  (Gould  et  al.,  2004), 
in  situ  ocean  observation  is  still  very  sparse  over  the  vast  ocean. 
Especially,  when  operational  modeling  and  forecasting  are  con¬ 
cerned,  it  is  often  difficult  to  And  in  situ  observations  near  the  ty¬ 
phoon  track  (Lin  et  al.,  2005,  2008,  2009b;  D’Asaro  et  al.,  2011; 
Pun  et  al.,  201 1 ).  Satellite-derived  OTS  thus  provides  critical  com¬ 


plementary  information  because  it  can  provide  synoptic  OTS  maps, 
instead  of  point-wise  in  situ  observations. 

To  obtain  accurate  forecast  of  typhoon  intensity  change,  as 
accurate  as  possible  during-typhoon  SST  cooling  and  enthalpy 
fluxes  are  needed,  because  ocean  is  the  energy  source  for  typhoon’s 
intensification  (Emanuel,  1999;  Bender  and  Ginis,  2000;  Emanuel 
et  al.,  2004).  Pre-typhoon  OTS  profiles  along  the  typhoon  track 
are  thus  needed  as  input  to  predict  the  during-typhoon  SST  cooling 
and  enthalpy  fluxes.  In  the  often  absence  of  in  situ  profiles,  satel¬ 
lite-derived  profiles  can  be  used.  Pig.  14  compares  the  performance 
of  the  two  satellite-derived  methods,  using  the  co-incident  and  co¬ 
located  Argo  profiles  depicted  in  Pig.  7  as  the  truth.  Using  these 
profiles  as  input  to  an  ocean  mixed  layer  model  (Mellor  and  Yam- 
ada,  1982),  the  performance  is  assessed.  Here,  the  drag  coefficient 
(Q)  used  in  the  SST-cooling  prediction  is  the  high  wind  coefficient 
(suitable  under  typhoon  condition)  from  Powell  et  al.  (2003).  All 
simulations  are  performed  under  typhoon  translation  speed  (i.e., 
Uh)  of  5  m  s“^  because  it  is  the  averaged  Uh  observed  in  the  WNPO 
during  typhoon  intensification  (Lin  et  al.,  2008,  2009b).  Typhoon 
wind  forcing  from  35-65  m  s“^  with  an  interval  of  5  m  s“^  is  used, 
so  that  all  possible  wind  speeds  from  Categories  1  to  5  conditions 
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Fig.  14.  (a)  SST  cooling  and  (b)  enthalpy  flux  comparisons  between  using  the  TLM  (red  curves),  REGWNP  (green  curves),  and  Argo  (black  curves)  profiles  as  input  to  the 
mixed-layer  model.  The  results  shown  in  (a)  and  (b)  are  the  average  of  16  profiles  (Fig.  7)  runs.  The  Saffir-Simpson  typhoon  scale  from  Categories  1-5  based  on  the  10-min 
maximum  sustained  winds  is  depicted.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


Table  2 

SST  cooling  assessment  from  16  profiles  (Fig.  7)  as  input  to  the  ocean  mixed  layer  model.  Unit  is  in  °C. 


Wind  speed  (m  s  ^) 

35 

40 

45 

50 

55 

60 

65 

Average 

TLM 

1.09 

1.44 

1.62 

1.72 

2.02 

2.30 

2.54 

1.82 

REGWNP 

1.47 

1.94 

2.18 

2.33 

2.75 

3.18 

3.58 

2.49 

Argo 

1.76 

2.30 

2.56 

2.71 

3.15 

3.58 

3.99 

2.86 

can  be  simulated.  Further  detail  in  the  methodology  can  be  found 
in  Lin  (2012). 

After  obtaining  the  SST  cooling,  the  enthalpy  flux,  i.e.,  the  sensi¬ 
ble  (Qs)  plus  latent  (Ql)  heat  fluxes,  under  typhoon  condition  can 
be  estimated  using  the  bulk  aerodynamic  formulas  (Clone  and  Uhl- 
horn,  2003;  Jacob  and  Shay,  2003;  Powell  et  al.,  2003;  Black  et  al., 
2007;  Lin  et  al.,  2009b,  2013)  as  follows. 


(ls  =  C„W{Ts-Ta)p,Cpa, 

(6) 

Ol  =  CEW{q,  -  qJp„L„a, 

(7) 

where  Ch  and  Ce  are  the  exchange  coefficients  of  sensible  and  latent 
heat,  respectively;  W  is  the  wind  speed;  Ts  and  Ta  are  SST  and  near¬ 
surface  air  temperature,  respectively.  Ts  is  the  during-typhoon  SST, 
calculated  as  initial  SST  minus  SST  cooling;  Ta  is  taken  as  a  constant 
of  27  °C;  Qs  and  Qa  are  surface  and  air  specific  humidity,  respec¬ 
tively;  and  pa,  Cpa,  and  Lva  are  air  density,  heat  capacity  of  the  air, 
and  latent  heat  of  vaporization,  respectively. 

Fig.  14  shows  that  for  all  wind  speeds,  the  REGWNP  (green 
curves)  predicts  more  accurate  during-typhoon  SST  cooling  and 
enthalpy  fluxes  than  the  TLM  (red  curves).  It  also  shows  that 
the  TLM  underestimates  the  during-typhoon  SST  cooling 
(Fig.  14a)  and  results  in  significant  over-estimation  in  the  en¬ 
thalpy  fluxes  (Fig.  14b).  Table  2  quantifies  the  performance.  It 
can  be  observed  that  on  average  the  during-typhoon  SST  differ¬ 
ence  between  the  REGWNP  and  Argo  is  limited  to  0.37  °C,  but 
the  average  difference  between  the  TLM  and  Argo  is  1 .04  °C. 


Table  3 

Enthalpy  flux  assessment  from  the  SST  cooling  of  16  profiles  (Fig.  7).  Unit  is  in  W  m~^. 


Wind  speed  (m  s  ^) 

35 

40 

45 

50 

55 

60 

65 

Average 

TLM 

479 

486 

501 

600 

604 

565 

492 

532 

REGWNP 

367 

338 

327 

380 

334 

220 

73 

291 

Argo 

348 

314 

301 

342 

269 

145 

-25 

242 

Since  the  average  SST  cooling  using  Argo  floats  is  2.86  °C,  the 
corresponding  error  for  the  REGWNP  would  be  13%  (0.37  °C/ 
2.86  °C)  and  36%  for  the  TLM  (1.04  °C/2.86  °C)  (Table  2).  As  for 
the  enthalpy  fluxes,  the  average  difference  between  REGWNP 
and  Argo  is  49  W  m“^,  but  the  difference  between  the  TLM 
and  Argo  increases  to  290  W  m“^.  Since  the  averaged  enthalpy 
flux  from  the  Argo  is  242  W  m“^,  the  equivalent  error  between 
REGWNP  and  Argo  would  be  20%  (49  W  m-^1242  W  m'^),  while 
the  error  between  the  TLM  and  Argo  is  120%  (290  Wm“^/ 
242Wm-2)  (Table  3). 

Exploring  the  possible  reason  for  the  poor  performance  of  the 
TLM  derived  profile,  one  possibility  is  the  coarse  vertical  resolu¬ 
tion  (Fig.  6).  Since  only  two  isotherms  are  used  to  characterize 
OTS;  but  for  SST  cooling  estimation,  a  more  complete  vertical 
ocean  thermal  structure  is  required.  Therefore,  interpolation 
and  extrapolation  are  involved.  In  other  words,  since  the  TLM 
only  provides  D20  and  D26,  the  rest  of  the  thermal  structure 
is  obtained  through  interpolation  and  extrapolation  of  D20  and 
D26.  As  shown  by  Fig.  15,  this  can  lead  to  significant  deviation 
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Fig.  15.  Comparisons  of  (a)  D26,  (b)  D24,  (c)  D20,  (d)  D18,  (e)  D16  and  (f)  D14  between  the  REGWNP  (green  dots)  and  TLM  (red  dots)  and  Argo  observations  in  September  of 
2010.  D24,  D18,  D16  and  D14  are  obtained  by  interpolation  or  extrapolation  for  the  TLM.  Sample  number,  RMSD,  correlation  and  bias  with  respect  to  Argo  are  shown.  (For 
interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


from  reality  (see  the  red  dots).  As  for  the  REGWNP  method,  since 
each  isotherm  is  derived  individually,  the  accuracy  is  much  high¬ 
er  and  closer  to  the  Argo  observations  (see  the  green  dots  in 
Fig.  15). 

7.  Conclusion 

Though  the  importance  of  ocean’s  subsurface  thermal  structure 
in  tropical  cyclone’s  intensity  change  has  been  known  since  1970s, 
it  was  difficult  to  obtain  subsurface  information  over  the  vast 
oceans  with  limited  in  situ  point-wise  observations.  Only  since 
the  launch  of  satellite  altimetry  in  1993  and  the  pioneering  work 
(i.e.,  the  TLM)  of  Goni  et  al.  (1996)  and  Shay  et  al.  (2000),  it  was 
possible  to  operationally  derive  synoptic  ocean  subsurface  thermal 
structure  using  altimetry  for  typhoons  intensity  forecasting  and 
research. 

However,  as  sciences  progress  it  is  desirable  to  improve  vertical 
characterization  of  OTS,  after  all  the  use  of  only  two  isotherms  in 
the  TLM  is  indeed  too  coarse.  With  the  deployment  of  Argo  floats 
in  2000,  there  has  been  a  quantum  jump  in  ocean  subsurface  sam¬ 
pling.  Therefore  it  has  become  possible  now  to  combine  this  large 
amount  of  in  situ  profiles  with  satellite  altimetry  to  regression  and 
develop  new  generation  of  derivation.  Based  on  this  method,  we 
derive  a  maximum  of  26  isotherms  (from  D4  to  D29)  to  improve 
the  vertical  characterization  (Fig.  6).  With  these  high  resolution 
profiles,  the  shape  of  ocean  subsurface  thermal  structure  can  be 
well-characterized.  Using  an  independent  set  of  more  than  7000 
Argo  profiles  for  validation,  it  is  found  that  this  new  generation 
of  derivation  is  of  good  accuracy  for  all  isotherms  with  error  mostly 
below  15%  (Fig.  9c).  Because  many  more  isotherms  can  be  used, 
evident  improvement  in  subsequent  parameters  derived  or  esti¬ 
mated  using  these  satellite-derived  profiles  is  achieved.  Improve¬ 
ment  is  found  in  all  related  parameters,  including  TlOO,  UOHC, 


during-typhoon  SST  cooling,  and  air-sea  enthalpy  fiux  estimation 
(Figs.  7-15). 

Further  validation  on  the  regional  dependence  found  that  the 
performance  is  better  over  the  southern  part  of  the  WNPO  (i.e., 
at  the  south  of  26°N),  co-located  with  the  main  development  re¬ 
gion  of  typhoons  in  the  WNPO.  Over  the  northern  region,  the  per¬ 
formance  degrades  by  about  20%  in  error  (i.e.,  30-40%).  Similar 
northern  degradation  in  performance  is  also  found  in  the  tradi¬ 
tional  TLM  but  is  much  more  severe.  As  in  Fig.  11,  the  error  found 
in  the  TLM  for  D20  can  be  more  than  100%.  Therefore,  over  the 
northern  WNPO  (i.e.,  at  the  north  of  26°N,  especially  over  the 
NEZ),  it  should  be  cautious  in  using  the  satellite-derived  OTS 
information. 

Finally,  as  this  new  generation  of  derivation  can  be  extended  to 
D4,  i.e.,  ~1000m,  the  application  of  satellite-derived  OTS  is  not 
limited  to  typhoon-ocean  interaction  application,  but  can  be  ex¬ 
tended  to  new  oceanographic  issues  because  information  from 
deeper  ocean  can  be  obtained. 
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